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The term 'tomography' is commonly applied to the idea of studying properties of a medium by 
the modifications this medium induces to a known probe propagating through it. In the context of 
ultrarelativistic heavy-ion collisions, rare high transverse momentum (pr) processes taking place 
alongside soft bulk-matter production can be viewed as a tomographic probe as long as the energy 
scales are such that the modification of high pr processes can be dominantly ascribed to inter- 
actions with the medium during the propagation of partons. Various high pj observables have 
been suggested for tomography, among them hard single hadron suppression, dihadron correla- 
tions and 7-hadron correlations. In this paper, we present a model study of a number of different 
observables within the same calculational framework to assess the sensitivity of the observables 
to different properties of the medium and discuss the prospects of obtaining tomographic infor- 
mation. 
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1. Introduction 



The expression 'jet tomography' is often used to describe the analysis of hard perturbative 
Quantum Chromodynamics (pQCD) processes taking place inside the soft matter created in an ul- 
trarelativistic heavy-ion collision. Such processes, which are well understood in p-p collisions, can 
be viewed as a known probe as they take place before any formation scale of a soft medium. Thus, 
only the subsequent propagation of partons through the soft medium (and possibly hadronization, 
although at sufficiently high pj the hadron formation length is larger than the medium extension) 
is sensitive to interactions with the medium. Hence, any modification of hard processes embedded 
in a medium potentially carries tomographic information about the medium properties. 

In particular the experimental focus is on the nuclear suppression of hard hadrons in A-A col- 
lisions compared with the scaled baseline from p-p collisions, which is expected due to interactions 
of a hard parton with the soft medium (see e.g. [[T]]). However, the nuclear suppression factor 

is a rather integral quantity, arising in model calculations from a convolution of the hard pQCD 
vacuum cross section do^c^^^^ for the production of a parton /, the energy loss probability 
Pf{AE) given the vertex position and path through the medium and the vacuum fragmentation 
function Dy?^^(z,/i|), as schematically expressed 

/ 

where 

dO^r^^'' = Y^fi/Aixi , e') ® /;M(^2, e') ® a,-y^/+^ (1-3) 
ijk 

Here, fi/Ai^M^) denotes the nuclear parton distribution function which depends on the parton 
momentum fraction x and the hard momentum scale and Ojj^f+k is the the partonic pQCD 
cross section. 

Eq. ( pT^ ) has to be properly averaged over all possible vertices distributed according to the 
nuclear overlap Taa and all possible paths through the medium. In [0] we have argued that one can 



factorize this spatial averaging from the momentum space formulation Eq. ( |1.2[ ) and thus define the 
geometry-averaged energy loss probability {P{AE,E))tj^^. Raa can thus be viewed as providing 
constraints for the form of {P{AE ,E))t^. 



2. Calculational framework 



Any model for medium modifications of a hard process must contain three major ingredients: 
The hard pQCD process, the bulk matter evolution for which we either use a hydrodynamic [||] or 
a parametrized evolution model [Q] and the energy loss probability distribution given a hard parton 
path through the soft medium [Q]. 

The primary hard process is calculated in leading order pQCD under the assumption that the 
transverse momentum scale is large enough so that hadronization takes place outside the medium 
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and that the produced leading hadron can be assumed to be collinear with its parent parton. The 
calculation, when supplemented by a K-factor, agrees well with hard hadron production measured 
in p-p collisions. In particular, the AKK set of fragmentation functions ^ also gives a satisfactory 
description of proton production whereas the older KKP set does not. Explicit expressions for 
the hard process calculation can be found e.g. in [||]. 

The interaction of the hard parton with the soft medium is calculated using the radiative energy 
loss formalism of ||5|]. If we call the angle between outgoing parton and the reaction plane 0, the 
path of a given parton through the medium (§(t) is specified by (ro,0) and we can compute the 
energy loss probability P{AE)path for this path. We do this by evaluating the line integrals 

poo OCX) 

COc(ro,^)= d^^m and {qL)M) = d^q{^) (2.1) 

JO JO 

along the path where we assume the relation 

<7((^) =is:-2-£^/'*(,^)(coshp-sinhpcosa) (2.2) 

between the local transport coefficient q{t,) (specifying the quenching power of the medium), the 
energy density e and the local flow rapidity p with angle a between flow and parton trajectory 
[^, [T^. Here cOc is the characteristic gluon frequency, setting the scale of the energy loss probability 
distribution, and {qL) is a measure of the path-length weighted by the local quenching power. We 
view the parameter ^ as a tool to account for the uncertainty in the selection of and possible 
non-perturbative effects increasing the quenching power of the medium (see discussion in [^]) and 
adjust it such that pionic Raa for central Au-Au collisions is described. Using the numerical results 
of [^], we obtain P{AE;(Oc,R)path for (Oc and R = 2(0^ / (qL) as a function of jet production vertex 
and the angle (/). 

The information about the soft medium is contained in the local energy density e(i§) and the 
flow rapidity p(i^). These parameters are obtained from dynamical evolution models which are 
tuned to describe a large body of bulk matter observables [^, Details of the evolution models 
including contour plots of their time evolution can be found in [^. In the following, we mainly 
illustrate three scenarios: A hydrodynamical evolution of matter ('Hydrodynamics'), the best fit 
to soft hadronic pj spectra and HBT correlation data of the parametrized evolution model ('Box 
density') and the hydrodynamical model under the assumption that only the partonic evolution 
phase leads to energy loss ('Black core'). Since in all models Raa for central Au-Au collisions is 
described by construction via a fit of K, the latter model implies that K takes large values and the 
evolution exhibits a very black interior region and a dilute hadronic halo which does not induce 
energy loss at all, quite different from the other models. 



3. Single Hadron Suppression 

Since Raa does not contain any spatial information, the production vertices of hard partons 
and their path through the medium have to be averaged out. Hard vertices {xo,yo) are distributed 
according to a probability density 
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Figure 1: Left panel: Trial energy loss distributions {P{AE))tj^ for various scenarios of jet energy loss in 
the medium (see text and for details). Right planel: Raa as calculated from the trial distributions shown 
on the left hand side. 



where b is the impact parameter. The thickness function is given by the nuclear density (r, z) as 
TA{r) = J dzPA{r,z)- Hence, given the energy loss probabihty distribution Pf{AE)path for a given 
path through the medium, we obtain 

{PfiAE,E))T^ = —J^^ Jdxo JdyoP{xo,yo)Pf{AE)path. (3.2) 

Before we proceed to calculate this quantity, let us illustrate the sensitivity of Raa to details of 
{Pf{AE))T^ (and hence the potential for tomographic information) by inserting trial distributions 
into the folding integral Eq. (pr2|). These trial distributions are shown in Fig. left panel, the 
resulting Raa is shown in the right panel and compared with the PHENIX data for pions [ jlTl ] (see 
also for details). 

It is apparent from the figure that despite strong differences in the functional form of {P{AE))t^ 
all distributions describe the measured Raa reasonably well above some minimum pj. The notable 
exception is the case of a constant fractional energy loss in which Raa drops as a function of pr, 
which does not seem to capture the overall trend well. It has to be concluded that Raa does not 
exhibit great tomographic capability beyond a single overall energy loss scale (the numerical value 
of which moreover is different for each model). This may explain why different calculations extract 
rather different quenching properties of the medium from fits to Raa- 

However, while the curves are reasonably similar over the kinematic range shown here, they 
do show differences in details which unfortunately cannot be resolved within the current data pre- 
cision. Thus, there is some reason to suspect that either increased ecperimental statistics or a larger 
accessible kinematic range may provide more stringent constraints for the energy loss distribution. 

In Fig. ^left panel we show the calculated Raa using the procedure outlined above to determine 
{P(AE))t^ instead of a trial ansatz. Once K is adjusted, the result does not exhibit strong sensitivity 
to the underlying medium evolution model (we do not show the results for all different medium 
evolutions here as the curves are difficult to distinguish), again confirming that Raa has very limited 
tomographic capability in the RHIC kinematic range. 

If the AKK fragmentations are used for computation for which the baseline process of 
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Figure 2: Left panel: Raa for pions and protons as compared to the measured {Rcp) data [ |1 1| , |12[ ] for 
RHIC conditions Right planel: Model predictions for Raa at the LHC based on two scenarios which describe 



the data at RHIC 114|. 



proton production in p-p collisions is roughly under control (in detail, AKK seems to overpredict 
the process by about a factor 2), the calculation of Raa for both pions and protons agrees well with 
the data [|^. This is not a trivial result, as in the AKK fragmentation scenario proton production 
is gluon-dominated whereas pion production is not, hence the difference between pion and proton 
production should reflect the different energy loss properties of quarks and gluons. In the calcula- 
tion as presented here, the rather small difference between proton and pion suppression is caused 
by the fact that gluon suppression is already in a saturated regime — increasing the quenching 
power of the medium further induces only a small change in the gluonic Raa 1 13]- 

In Fig. ^ right panel we show the extrapolation of the pj dependence of Raa to LHC energies 
based on the hydrodynamical scenarios which describe the data at RHIC [^,14]. While there 
is some uncertainty associated with the extrapolation of the nuclear parton distribution function 
(NPDF [||] vs. EKS98 this is a small effect, and it becomes indeed apparent that with 

the extended kinematic lever-arm of LHC the different properties of the two scenarios (dense core 
and dilute halo vs. more evenly distributed quenching power) can clearly be distinguished. The 
results here differ from a previous calculation presented in [p^. The improvement of the present 
calculation over the previous work is chiefly in the use of a dynamically evolving soft medium 
instead of a static cylinder ansatz and in using Eq. (3J) for the primary vertex distribution as 
compared to a homogeneous distribution. 

Let us illustrate the differences induced by the spatial distribution of the quenching power by 
studying the geometry of single hadron suppression directly in the model. In Fig. ^ we show the 
probability density of finding the primary pQCD vertex leading to an observed hadron above 8 GeV 
in pt- It is evident (and quite expected) that emission occurs predominantly close to the near side 
surface of the medium. However, the degree to which surface emission is realized is quite different 
in all three models. Clearly, the strong suppression from the core region of the black core scenario 
repels the distribution much more from the center than the more even distribution of the other two 
scenarios. Thus, surface emission is not a property of a particular energy loss formalism, but arises 
from the interplay of energy loss formalism with the underlying geometry and evolution of the soft 
medium. 
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Figure 3: Probability density for finding a hard vertex in the transverse plane in 200 AGeV Au-Au colli- 
sions leading to an observed hadron above 8 GeV transverse momentum, shown for three different medium 
evolution scenarios (see text). In all cases, the hard hadron propagation defines the —x direction. All contour 
intervals are linear 



4. Dihadron suppression 

We can make use of the sensitivity of the vertex distribution of single hadron suppression to 
the medium evolution by considering dihadron suppression. In a back-to-back event, the second 



hadron propagation path is not averaged over the initial overlap Eq. (3.1) but over a conditional 
probability distribution given a valid trigger, i.e. over the distribution shown in Fig. ^ (which 
is quite different from the overlap). Thus, even if two model evolutions lead to identical Raa, 
this does not mean that they would produce the same dihadron correlation pattern. We call this 
conditional probability distribution given a high-pj- near-side trigger in the following {P{AE))jr 
and investigate its capability to obtain tomographical information. 

For computational purposes, we employ a Monte-Carlo (MC) simulation of the experimen- 
tal trigger condition, followed by the simulation of the away-side parton intrinsic-^j- smearing, 
propagation, energy loss and fragmentation. The procedure is described in detail in 

In Fig. ^ we compare the yield per trigger on the near and away side for different medium 



models with the data obtained by the STAR collaboration |19p. Within errors, the near side 
yield per trigger is described by all the models well. There is no significant disagreement among 
the models. The model calculations appear significantly more different if we consider the away 
side yield. Here, results for the 4-6 GeV momentum bin differ by almost a factor two. However, 
none of the model calculations describes the data in this bin. This is in fact not at all surprising 
as below 5 GeV the inclusive single hadron transverse momentum spectra are not dominated by 
pQCD fragmentation and energy losses but, rather, by hydrodynamics possibly supplemented with 



recombination [gO|, |21|] type phenomena. For this reason, the ratio Raa at /j^ < 5 GeV cannot be 
expected to be described by pQCD fragmentation and energy losses, either. 

This is clearly unfortunate, as the model results are considerable closer to the experimental 
result in the 6+ momentum bin on the away side and hence our ability to discriminate between 
different models is reduced. Since at this large transverse momenta the pQCD fragmentation -i- 
energy losses dominate the singe hadron spectrum, we expect that the model is able to give a valid 
description of the relevant physics in this bin: Not only is Raa well described by the data, but also 
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Figure 4: Yield per trigger on the near side (left panel) and away side (right panel) of hadrons in the 4-6 Ge V 
and 6+ GeV momentum bin associated with a trigger in the range 8 GeV < pj < 15 GeV for the different 
models of spacetime evolution as compared with the STAR data JT^ |l9|. The individual data points have 
been spread artificially along the x axis for clarity of presentation. 
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Figure 5: Yield per trigger on the near side (left panel) and away side (right panel) of hadrons in the 4-6 
GeV and 6-1- GeV momentum bin associated with a trigger in the range 12 GeV < pr < 20 GeV for the 
different models of spacetime evolution. The individual data points have been spread artificially along the x 
axis for clarity of presentation. 



the contribution of recombination processes to the yield is expected to be small [20|. Thus, as it 
stands, only the black core scenario can be ruled out by the data, the box density with Bjorken 
expansion seems strongly disfavoured but still marginally acceptable. 

Thus, as it stands, the kinematic window to study dihadron correlations in a perturbatively cal- 
culable region is not enough to exploit the difference between {P{AE))Tr and {P{AE))}i^ and thus 
to obtain detailed tomographic information. However, the situation may improve for an increased 
kinematical window in the region where pQCD -i- fragmentation can be applied. In order to test 
this, we redo the MC simulation with trigger hadrons in the range between 12 and 20 GeV. 

The distribution after fragmentation into hadrons in bins of 2 GeV width in the perturbative 
region is shown in Fig. ^ for the near side (left panel) and away side (right panel). It is again appar- 
ent that within errors all models agree in the expected near side yield. The momentum spectrum 
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Figure 6: Probability density for finding a hard vertex in the transverse plane in 200 AGeV Au-Au collisions 
leading to an event with both an observed near side hadron above pT — 8 GeV (defining the — x direction) 
and an away side hadron with pj > A- GeV for three different medium evolutions (see text). All contour 
intervals are hnear. 
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Figure 7: Left panel: Probability density for finding a hard vertex in the transverse plane in 5.5 ATeV Pb-Pb 
collisions leading to a near side hadron with pj > 25 GeV propagating into the — x direction for the LHC 
hydrodynamical model prediction. Right: Probability density requiring in addition an associated away side 
hadron with pr > 10 GeV momentum. All contour intervals are linear. 

of the away side exhibits considerably more structure. Several of the scenarios can now be clearly 
told apart in bins in the perturbative region. For example the Ta and the box density (which have 
virtually identical (/'(AE))^^) show almost a factor two difference in the 10-12 GeV momentum 
bin. As we have seen above in the case of the LHC extrapolation, is evident again from the analysis 
that having a larger lever-arm in momentum is needed to get access to tomographic information. 

Finally, let us discuss the geometry of dihadron suppression. In Fig. ^ we show the probability 
density of finding a hard vertex leading to a high pr near side trigger and a correlated associated 
hard away side hadron. Here, clear differences between tangential emission in the case of a dense 
core and production across the whole volume become apparent. 

In Fig. ^ we also show the geometry of single hadron and dihadron suppression for LHC 
conditions in central Pb-Pb collisions. For a 25 GeV trigger hadron, we expect some degree of 
surface emission (note that the dihadron production distribution is somewhat repelled from the 
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center) but no strong tangential emission. 
5. A simple model 

As we have seen, quite a general class of models predict a small rise of Raa with pj at RHIC 
and a more pronouned one at LHC. Assuming RHIC kinematics, Raa is rather insensitive to de- 
tails of the energy loss probability distributions, at LHC the sensitivity is considerably enhanced. 
Likewise, dihadron correlations become more sensitive to the medium density distribution if the 
kinematic range is increased. In the following, let us try to illustrate that all these observations can 
be understood from simple considerations. 

Quite generally, energy loss probability distributions can be decomposed as 



where T is a transmission term describing a parton penetrating through the medium without energy 
loss, 5 is a shift term which characterizes partons emerging from the medium after a finite energy 
loss and A is an absorption term describing partons which have been shifted in energy so much that 
they become part of the soft medium. 

Let us now assume a power law for the parton spectrum at RHIC and LHC ^ \/pj with 
nRHic > niHC- Energy loss AE then changes this spectrum to l/(pr + AE)", thus Raa in this 
simple model can be obtained from 



It is evident from the expression that Raa at given pT is equal to the transmission term T plus 
a contribution which is proportional to the integral of {P(AE))t^ from zero up to the energy scale 
E,„ax of the parton, seen through the filter of the steeply falling spectrum. Raa grows with pj since 
Eynax grows linearly with pj. However, at RHIC conditions the characteristic scale <x>c of the energy 
loss probability distribution is far above Emax^ thus the growth is slow and Raa is dominated by T, 
rendering it almost a constant. Since tomographic information is mainly contained in the shift term 

5, the apparent insensitivity of Raa to assumptions about the medium can be understood. 

This is very different at LHC where E^ax ~ o>c (since Emax grows linear with pj but (Oc grows 
with the entropy density and hence much slower) and a pronounced contribution of the shift term 
can be probed. Here, a rise of Raa with pj is expected, along with a greater tomographic sensitivity. 

6. Conclusions 

We have investigated the capability of single and dihadron suppression to provide tomographic 
information about the soft medium created in ultrarelativistic heavy-ion collisions. We have argued 
that at RHIC kinematics, the nuclear suppresssion factor Raa is not very sensitive to the medium 
evolution. While dihadron suppression, due to its different geometrical averaging, exhibits in prin- 
ciple more sensitivity to medium properties, unfortunately the present data situation allows only to 
rule out a very pronounced difference between a strongly absorbing core and a dilute halo. This 



{P{^E))T^^Tr) =T5{AE) + S- P{AE) +A-5{AE-E) 



(5.1) 
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insensitivity can be traced back to the fact that (Oc, the intrinsic scale for energy loss is much higher 
than EfficLx, the accessible parton energy at RHIC. 

However, when going to LHC energies, this condition no longer holds. Raa becomes dom- 
inated by partons being shifted in energy, and tomographic information can be recovered even 
from the pT dependence of Raa- Dihadron correlations and other measurements, such as 7-hadron 
correlations which provide a monochromatic source of hard quarks in the medium or Raa vs. 
reaction plane [p2|], which allows for a systematic variation of in-medium pathlength, may pro- 
vide additional information such that a multi-pronged approach to jet tomography finally becomes 
feasible. 
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